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USER'S GUIDE AND DESCRIPTION OF THE PERTURBATION ` 
CODE MODULES DGRAD and TPERT 


Е. T. Tomlinson and В. A. Lillie 


ABSTRACT 


Two code modules, TPERT and DGRAD,are described, which 
when used in conjunction with the DOT IV discrete ordinate 
neutron transport code provide the user with the capability 
of performing perturbation calculations. The modules and 
their input are described. А sample problem is presented 
to aid in executing these modules. 


1. INTRODUCTION 


The movement and/or replacement of one material by another in a 
reactor core will cause a change in the system multiplication factor, 
keff: The magnitude of this change can be calculated either directly 
by recalculation of keff for the perturbed system or through the use of 
perturbation techniques. The methods incorporated in the code modules 
described in this report include first-order and exact perturbation 
theory based on both transport and diffusion theory approximations to 
the neutron transport equation. These modules were written as ancillary 


codes to the DOT IV discrete ordinates transport code. All input to 


these models is done uSing the free-form input routine FIDO. 
2. TPERT 


A perturbation code, TPERT, has been written to obtain exact and 
first-order reactivity changes. The purpose of this section is to 


provide a brief review of the expressions evaluated by TPERT and to 
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provide a detailed description of the required input. ТРЕКТ couples with 
VIP! which generates adjoint-forward flux tables, i. e., ф*ф tables, 
using DOT-III? scalar flux tape information. (A slight modification of 
VIP will allow the use of DOT-IV scalar flux tapes.)? The general scheme 
consists of passing the spatially integrated ф*ф tables from VIP into 
TPERT which in turn folds the 4*6 tables with the pertinent cross section 
changes. 

The reactivity associated with a change in reactor composition, i.e., 


a perturbation, may be obtained from: 


Ap = <ф*,(АЕ/К-АЕ+А5$)ф7>/<ф*,Е^ф7> (1) 


where 
Ар = reactivity change 
AF = change in fission operator 
AL = change in total cross section operator 


AS = change in scattering operator 


К = unperturbed state eigenvalue 
„* = unperturbed adjoint flux 
° = perturbed forward flux 
F“ = F + AF = perturbed fission operator. 


The brackets, i.e., <>, in Equation (1) indicate integration over 
space, energy and angle. Equation (1) is normally called the exact 
perturbation equation since there are no approximations employed in its 
development. The first-order, or linear, approximation to Equation (1) 
may be obtained by simply substituting the unperturbed forward fiux 9 


and the unperturbed fission operator Р for 4^ and F^, respectively. 


This yields: 
Ар = , (AF/k-Aztas)e»/«o*,Fé» . (2) 


The evaluation of the integrals in Equations (1) and (2) is carried 
out in two steps. The first step is performed by VIP and results in the 
generation of (141), where L is the order of scattering anisotropy, $*4 
tables for each region in which material changes occur or in which 
fissionable nuclides appear. These tables are ordered in a manner similar 
to ANISN fcrmatted cross-section tables. 

The terms which fold with the total cross section in TPERT are 


located in position IH7. These terms have the form“: 


241 : : 
(-1) (253138, 5, j 9 -— Vi (3) 


where 


Jo ming = the DOT 2m" adjoint flux moment in space cell i 


and group g 
Jimi g = the DOT forward flux moment. 
Vi - volume of jth space cell. 


The terms which fold with the scattering and fission operators are 
stored in positions IHS through IHS+G-] where G represents the total 
number of groups. These terms have the form": 

j* 


£ ,m,1,g ја mg Ч. (4) 


For a given group g, the subscript g^ may take on two vaiues. In po- 


sitions IHS through IHS*g-1: 

g^ = 9-К+1 for k-1,...,9g 
and in positions IHS+g through 1Н546-1: 

97 = k for К=9+1...,6 . 


The second step in the evaluation of the integrals in Equations (1) 
and (2) is carried out by TPERT. This step consists of folding the terms 
represented by Equations (3) and (4) with the appropriate cross-section 


changes for each region. The folding process is described by the 


following: 
АЕ E Egg 
. K Òr >}, Xg Kk 10, 9,9 ° (54) 
= * = 
Oe, areh Jorg т, (55) 
апа 
, ASO). = AL " 5 5b 
< Pr LLL 82,9 +g 2,9,9 (90) 
where 
4g = fraction of fission neutrons born in group g, 
ду гдг = change in nu times the fission cross section in 


group g 
Ас = change ín total cross section ín group 9, 


А = спапде 1г pth order scattering coefficient for 


HN 
87.9 +9 
scattering from group g^ to group д. 


In Equations (5), 


material regions. 


the subscript r has been introduced to denote different 


The perturbation denominators, i.e., , F p$” and 


<ф*,Еф> are evaluated by simply substituting 6 or ур, into 


g 


Equation (5a) for ur Summing the region contributions to Ap com- 


pletes the reactivity calculation. 


TPERT INPUT 


Title Card (18M) 


1$$ IHT- 
185 - 


Са 
+ 
= 
e 
u 


ITOD 


Position of total cross section 

Position of self-scatter cross section 
Cross-section table length 

Number of energy groups 

Number of "VIP" regions (see Note 1) 

Order of scattering anisotropy 

Number of P, materials from tape (see Note 2) 
Logical unit number of VIP 4*4 tables 

Logical unit number of cross-section tape 


Direct access logical unit number for cross-section 
storage (see Note 3) 


Logical unit for scratch disk (diffusion only) 
0 - use perturbed VEg dencminator 


1 


use unperturbed У denominator 


0 - transport calculation (see Note 8) 


- 
1 


diffusion calculation 


Terminate data block 1 with a T. 
2** CHI (IGM) - Fission spectrum 


3** VELS (IGM) - Group inverse velocities - fill with 1.0's if 
inverse velocity folding is not desired (see Note 7) 


4$$ MATID (NMAT)- Cross-section tape ID's for Ру materials. 
Arrays 555 through 1055 apply to the following relationships for a given 
VIP region: 
£ unperturbed = Ul*z,4U2*;, 
= perturbed = P1*z,+P2*55 
555 MATI (НВ) - Ру material numbers on direct access storage for 
Үү cross sections by VIP region - use negative 
number if region has fissionable material but no 
perturbation. 
£** U) (NR) - Unperturbed multiplier for £) by VIP region 
7** pi (№)  - Perturbed multiplier for £) Бу VIP region 
855 MAT? (NR) - Р, material numbers on direct access storage for 
rj cross sections by VIP region - use a zero if 
only Li occurs in a region 


9** U2 (NR) - Unperturbed multiplier for z, by VIP region 
10** P2 (NR) ; Perturbed multiplier for z, by VIP region 


11** EIG - Unperturbed Koff 
BETA - Total delayed neutron fraction 
DNORM - Perturbation integral denominator - enter 0.0 


if code is to calculate 
Terminate data block 2 with a T. 
Note 1: Each region which contains a perturbation, i.e., change in 
cross section, or which contains fissionable material is a 
VIP region. 


Note 2: ТРЕКТ, as presently programmed, will only accept cross sections 
їп ANISN library format - either on tape or disk. 


Note 3: Direct access is used te store the cross sections in the same 
order as the ID's in the 455 array, e.g., if the ID of material 
M corresponds to the Nth entry in the 455 array, the Р, rumber 
for this material is N in the 555 and 8$$ arrays. 

Note 4: The number of words needed in COMMON/DATA in TPERT main is 
given by: 

(3*2* ITL) *16М+ММАТ * ( ISCT*2) * 10*NR*40 
Note 5: The number of tracks needed on logica; unit КОСЕ is given by: 
( (ITL*TGM-1)/1800*1 ) *RMAT* C ISCT*1) 


Note 6: To run multiple cases, simply supply a new title card and change 
the appropriate arrays. Be sure to include a terminator, e.g., 
T for both data blocks. 


Note 7: TPERT will evaluate the quantity «;*, E +> which is useful in 
certain time-dependent reactivity cal chlat ions. 


Note 8: То run diffusion cases, it is necessary to set ISC I-] in VIP 
ana input the flux iapes generated by DCRAD. 


3. DGRAD 


A module, DGRAD, was written to calculate the directional flux 
gradients from DOT III? diffusion theory flux tapes. 00Т-1У° scalar 
flux tapes can also be esed if they are preprc-essed in the RTFLUM 
module of the 0055 system. The module, DGRAD, in conjuncrion witn VIP’ 
and TPERT will allow the use of diffusion theory fluxes to obtain exact 
and first-order perturbation reactivity changes. 

In order to calculate the reactivity associated with changes in 
reactor compositions using diffusion theory, it is necessary to fold not 
only the scalar fluxes with the appropriate cross sections, but also the 


average flux gradients with the diffusion coefficients. Since DOT 


diffusion thecry'does not directly calculate these gradients, it was 


necessary to calculate the needed quantities external to the DOT code. 


Consider the typical mesn cell with ceil-centered fluxes 


e ч 1 
The flux on the left boundary of cell i, 3, $gi can be calculated if 
the current across this boundary is assumed to be continuous. 


p — ói. . $i 2 
220: BL "i-li вор 23-28 
1-1.) bey 1,3 5 


ипеге 0, 1 is the ditfusion coefficient for cell i,j. Solving for 


5 


8L 


„ 9111-17401 
BL Di 3j 41 + Dij А-1 


Similarly for the right boundary 


Ding ij 1 * 76 1,1 ^in 


“BR Ü; j ia В, Ч 


The R direction gradient can then be obtained as follows 


Similar arguments are used for the calculation of the Z direction 
gradients. 

On the outer boundaries of the problem the following conditions 
are imposed 


For reflected boundaries at 


EE 27903 jl 0 1 
90, 413 51207515] 
00,5 ^ P1, Di 0 7 Pi 

Мом, = PIM,j jM 97 omer 2 °, 
IN, ^ *IM,j ĉi % , 0% 
DIM+1,j bin. Di 1,0% 


For non- return boundaries the Flux is assumed to equal zero at ihe 


extrapolated boundary 


цагг 5 = 0.0 jl o; 9 = 0.0 
àg.; = 1-42089/түү 0 1.42089/2a 
90, ^ 01,3 9% * Di | 
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i-IM I. ,j = 0.0 13-48 : ЈМ = 0.0 
qug = 1-42089/: тр ИЕ ут 
Drwe1,j ^ 01м, ј О мн 7 D; JM 
DGRAD INPUT 


Title Card (20A4) 


1$$ IFLX 


0/1 Forward/Ad joint 

IM - No. of Radial Intervals 

JM  - No. of Axial Intervals 

MZ  - No. of Material Zones 

ITL - Cross-Section Table ! ength 

[НТ - Position of the Transport Cross Section 

МТ - Total Number of Materials 

IGM - No. of Energy Groups 
Left Boundary Condition 0=УАС 1=ВЕР 
- Right Boundary Condition 

IBT - Top Boundary Condition 


IBB Bottom Boundary Condition 

IP] 0/1 Diffusion/Hybrid Output (see Note 2) 
NIN Flux Input Unit 

NOUT - Flux Output Unit 

NSEC - Cross-Section Input Unit 
ISCRT - Scratch Unit 

IPRT - 0/1 Print/No Print Flux and Gradients 


Data Block must be terminated with a T. 


q** 
5** 


(IM+1) Radial Boundaries 
(JM+1) Axial Boundaries 


6$$ (IM*JM) Zone Number By Interval 


7$$ 


(MZ) Material By Zone 


Data Block must be terminated with a T. 


Note 1. 


Note 2. 


To run multiple simply supply a new title card and input all 
the arrays and T's. 


If 1Р1 is set equal to zero the scalar flux and the flux 
gradients are output for the TPERT diffusion perturbation 
calculation. If IP! is set equal to 1 the scalar flux and 
the partial currents are output for the VIP transport íntegral 
calculation. VIP is then run as if it were doíng a normal P, 
transport problem. The latter option allows for including 

the P, scattering effects and implicit calculation of neutron 
leakage worths. 
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